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Abstract. The three species ABC model of driven particles on a ring is generalized 
to include vacancies and particle-nonconserving processes. The model exhibits phase 
separation at high densities. For equal average densities of the three species, it is 
shown that although the dynamics is local, it obeys detailed balance with respect 
to a Hamiltonian with long-range interactions, yielding a nonadditive free energy. 
The phase diagrams of the conserving and nonconserving models, corresponding 
to the canonical and grand-canonical ensembles, respectively, are calculated in the 
thermodynamic limit. Both models exhibit a transition from a homogeneous to a 
phase-separated state, although the phase diagrams are shown to differ from each 
other. This conforms with the expected inequivalence of ensembles in equilibrium 
systems with long-range interactions. These results are based on a stability analysis 
of the homogeneous phase and exact solution of the hydrodynamic equations of the 
models. They are supported by Monte-Carlo simulations. This study may serve as 
a useful starting point for analyzing the phase diagram for unequal densities, where 
detailed balance is not satisfied and thus a Hamiltonian cannot be defined. 
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1. Introduction 

Systems driven out of equilibrium by an external field, such as a temperature or a 
pressure gradient or by an electric field, have been a subject of extensive studies in 
recent years (see [l], 0| for reviews of this topic). Often, such systems reach a current 
carrying steady state where detailed balance is not satisfied. The nature of the steady 
state typically depends on the details of the microscopic dynamics. This is in contrast 
to systems in thermal equilibrium where the equilibrium distribution is determined by 
the Hamiltonian, and is independent of the microscopic dynamics. Thus, for instance, 
in the Ising model, the same equilibrium state is reached in the thermodynamic limit 
irrespective of whether the dynamics is magnetization conserving (Kawasaki dynamics) 
or magnetization nonconserving (Glauber dynamics). 

Insight into the properties of nonequilibrium steady states has been obtained by 
detailed studies of simple models of particles on a lattice, where the particles are driven 
by either boundary terms or in the bulk. The behavior of these driven systems has been 
shown to be nontrivial, as they exhibit phenomena which do not occur in equilibrium. 
In particular, nonequilibrium steady states of driven systems have been found to exhibit 
long-range correlations, even for systems with strictly local dynamics (sUsl]. While 
attempting to understand the nature of these long-range correlations, several studies 
revealed similarities between some properties of driven nonequilibrium systems and 
equilibrium systems with long-range interactions. One example is the presence of phase 
transitions in one-dimensional driven models 

0-E3, which are known to exist also in 



equilibrium systems with long-range interactions. An interesting question, which is 
addressed in the present paper, is to what extent can studies of equilibrium long-range 
systems provide some understanding of properties of driven nonequilibrium systems. 

Equilibrium systems with long-range interactions are those where the two body 
potential decays at large distance, r, as l/r d+a , with —d < a < in d dimensions. In 
this case the total energy of a system of homogeneously distributed particles scales as 
yi-a/d w ftYi its volume V. Therefore, if — d < o < 0, the energy scales super-linearly 
with V, and the system is nonadditive. As a result of the nonadditivity the various 



statistical mechanical ensembles may become inequivalent [I2|-|20j. For example, it is 
well known that the specific heat of long-range interacting systems can be negative in 
the microcanonical ensemble. On the other hand within the canonical ensemble the 
specific heat is proportional to the energy fluctuations, Cy oc (E 2 ) — (E) 2 , and always 
positive. Similar effects can be found when comparing the canonical and grand-canonical 



ensembles of long-range interacting systems 21|. It has been suggested that differences 



between ensembles in these systems are usually manifested in the vicinity of first order 



phase transitions 19l. |22|-|24|. For recent reviews of systems with long-range interactions 



see e.g. 



25|-|29|. Finding similar features in nonequilibrium driven systems would thus 
provide a further link between the two classes of models. 

A particularly interesting model within which the interplay between drive and long- 
range interaction may be conveniently analyzed is the ABC model. This is a one- 
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dimensional driven model on a ring where three species of hard-core particles evolve 
under particle conserving local stochastic dynamics. The model has been introduced by 
Evans et al. 



30L l31| who observed that although the dynamics of the model is local, it 
exhibits long-range order characterized by phase separation of the three species. It has 
been found that in the special case where the average densities of the three species 
are equal, the steady state distribution of the ABC model obeys detailed balance 



30, 



31 



This 



with respect to an effective Hamiltonian with long-range interactions 
observation suggests that although detailed balance is not satisfied for unequal densities 
and the steady state cannot be expressed in terms of an effective Hamiltonian, some 
characteristic features of the long-range interactions may still be observed in the driven 
regime where the densities are unequal. The ABC model has been considered by Clincy 
et al. j32| in the weak drive regime, where the driving force scales as the inverse of the 
system size. In this limit the model has been found to exhibit a second order phase 
transition from the homogeneous state to the phase-separated state at some value of 
the drive. 

In this paper, we introduce a generalization of the ABC model allowing for vacancies 
and processes which do not conserve the number of particles. We demonstrate that even 
in the presence of vacancies and nonconserving processes, the model possesses detailed 
balance when the average densities of the three species are equal. In this case, too, the 
equilibrium steady state may be expressed in terms of an effective Hamiltonian with 
long-range interactions. This enables us to compare the steady states of the conserving 
and the nonconserving dynamics by studying the free energy of the two ensembles. 
Under the equilibrium conditions where the three densities are equal, the conserving 
dynamics lead to a steady state corresponding to the canonical ensemble, while the 
nonconserving dynamics results in a grand-canonical equilibrium state. Since in both 
cases the effective Hamiltonian is long-ranged, we find that the two ensembles yield 

33]. 



different phase diagrams. A brief account of this study is given in 

The results presented in this paper may be used as a starting point for studying 
the nonequilibrium regime, where the densities are unequal and the steady state cannot 
be expressed in terms of an effective Hamiltonian. Such studies of the original, particle 
conserving, ABC model in the case of unequal densities have shown that the second 
order transition found for equal densities, persists for unequal densities, and its position 
varies continuously with the densities 0, 32|. 



The paper is organized as follows: The ABC model is defined and its properties 
are summarized in Section [2j In Section [31 the ABC model is generalized to include 
vacancies and nonconserving dynamics. The phase diagrams for the conserving and 
nonconserving models under effective equilibrium conditions are derived in Sections H] 
and[5l respectively. In Section[6]the model is further modified to explore some nongeneric 
features, found in the nonconserving phase diagram. We end with some concluding 
remarks in Section [71 



4 



2. The ABC model 

The ABC model is a p rototypical model of driven systems exhibiting long-range order 
in one dimension 30|, |3l(. It belongs to a wider class of models exhibiting similar 



34, 35fl in the context of sedimentation 



features, such as the one studied by Lahiri et al. 
processes. The ABC model is composed of three species of particles, labeled A, B and 
C, occupying a one- dimensional periodic lattice of length L. Each site in the lattice is 
occupied by a single particle. The number of particles of each type is given by Na, Nb 
and Nc respectively, where Na + Nb + Nq = L. The model evolves under random 
sequential dynamics that conserves the particle numbers of the three species, whereby 
a site is chosen at random and is exchanged with its neighbor in a clockwise direction 
according to the following rates: 

AB^BA, BchcB, CA^AC. (1) 
ill 

For q = 1 the dynamics is symmetric, and thus obeys detailed balance. In this 
case the model relaxes to a homogeneous equilibrium state, in which all particles 
are evenly distributed throughout the lattice. On the other hand for g / 1, the 
model relaxes to a nonequilibrium steady state in which the particles phase separate 
into three distinct domains. For q < 1 the domains are arranged clockwise in the 
order AA . . . ABB . . . BCC . . . C, and counterclockwise for q > 1. This is a strongly 
phase-separated state in the sense that fluctuations result in broadening of the domain 
boundaries to a finite width, leaving the bulk of the three domains unmixed in the 
thermodynamic limit. Since the steady states corresponding to q > 1 and q < 1 are 
simply related by space inversion symmetry we take q < 1 throughout the paper. 

As a result of the dynamical asymmetry of the model, the steady state of a finite 
system generally exhibits nonzero currents of particles. The net current of, say, the A 
particles is determined by the difference between the rate at which an A particle diffuses 
to the right through the B domain, ~ q Ns , and the rate at which it diffuses to the left 
through the C domain, ~ q N ° , so that 

J A ~q NB -q Nc . (2) 

All currents vanish in the thermodynamic limit where L — > oo with Na/ L, Nb/L, Nc/ L 
kept fixed. Equation ([2]) implies that in the special case of equal average densities, 
Na = Nb = Nc, the steady-state currents vanish even for finite L. Thus, although the 
system is driven by asymmetric forces, it reaches a steady state which seems to have no 
irreversible currents of particles. 

This result suggests that at equal densities, detailed balance may be satisfied, so 
that the model evolves into an equilibrium steady state. Indeed, it has been shown that 



the model obeys detailed balance with respect to an effective Hamiltonian 30, 131J . This 
Hamiltonian possesses long-range interactions, despite the local nature of its dynamics 
(CQ). It is defined in terms of the microscopic configurations of the model, which consist 
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of the set {X{\ = {A iy B iy Cj}, % = 1, . . . , L, where 

I 1 if site i is occupied by an X particle , , 

I otherwise. 

The Hamiltonian is then given by 

L-l L-i 

ft ({X,}) = ^ ^ (AC l+j + B t A i+J + CiB i+j ) . (4) 

i=l j=l 

The interaction between particles in 7-L is long ranged, mean-field like. The total energy 
of the model scales super-linearly with its length {% ~ L 2 ) which is characteristic of 
systems with long-range interactions. The periodic boundary conditions imply that the 
model is translationally invariant. In Eq. (Pflh however, site 1 is arbitrarily chosen. One 
can check that the Hamiltonian indeed yields the same energy regardless of this choice, 
as long as the densities of the three species are equal. A manifestly translationally 
invariant form is obtained by averaging over all possible choices of site 1, leading to 

L L-l , 

u ax,}) = J2J2j (AjBi+fc + BiCi+k + CiAi+k] ' (5) 

t=l k=l 

with the periodic boundary condition, Xl+i = X{. The two representations were 
shown to yield the same energy for every microscopic configuration [3l|. In these 
representations the energy of the ground state, where the three species are fully 
separated, is L 2 /9. 

Using this Hamiltonian, the steady-state distribution of the ABC model with 
Na = Nb = Nc = L/3 is given by: 

P{{X i }) = ±q H «™ (6) 

where Zl = Yl{x } q n ^ Xi ^ is the partition sum. The fact that detailed balance is 
satisfied with respect to the Hamiltonian (j4j) can be verified by considering an exchange 
of two particles, say AB — > BA. According to Eq. (j3J), the resulting change in "H 
due to this exchange is +1, whereas the reversed process changes l-i by —1. Indeed, 
the exchange of two neighboring particles of any different species leads to an energy 
increment AH = ±1. This, together with the expression ([6]) for the distribution function 
and the transition rates (OQ) leads to detailed balance. 

The model has been shown to exhibit phase separation for any q ^ 1, while relaxing 
to a homogeneous state for q — 1. In order to study the phase transition between the 
two types of equilibrium states, Clincy et al. 32| considered the model with an In- 
dependent q . Taking the limit of weak asymmetry, where q — > 1 as L — > oo, the model 
was found to relax to one of the phases depending on the rate at which q approaches 1 
at large L. It has been shown that for equal densities and a transition rate of the form 
q = exp (—(3/L), the steady state in the thermodynamic limit is homogeneous for < (3, 
and inhomogeneous for f3 > f3 c with a second order transition at f3 c = 2tt\/3. Since 
P({Xj}) ~ e -P n ({ x i})/ L ) the parameter (3 can be regarded as the inverse temperature 
of the model and 1/(3 C as the critical temperature. 
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The ABC model has also been studied on an open interval by Ayyer et al. j36|. In 
this case the model exhibits detailed balance for arbitrary average densities of the three 
species, and the phase diagram of the model in the entire space of densities has been 
derived in the weak asymmetry limit. The mean field approximation has been shown 



to be exact in the thermodynamic limit 



il, l36l-l38| , and an ana 



ytic expression for the 



36 



density profiles in the phase-separated state has been obtained 

In the following sections we generalize the ABC model on a ring to include 
nonconserving processes and analyze the resulting phase diagrams in the weak 
asymmetry limit. 

3. Generalized ABC model: vacancies and nonconserving processes 

We now introduce a generalization of the ABC model, allowing a comparison of two 
alternative dynamics: particle-conserving and particle-nonconserving. We begin by 
introducing vacancies (labeled 0) into the lattice. Thus, each site may be occupied 
by a particle of either of the species A, B or C or may remain vacant, 0, with 
Na + Nb + Nq = N < L. Vacant sites are dynamically neutral, so that a particle of 
any species may hop to the left or to the right into a vacant site with equal probability. 
Hence, the following rule is added to the exchange rules in Eq. ([I]): 

X0 +± OX, (7) 
1 

where X — A,B, C. 

We proceed by introducing a nonconserving process, whereby triplets of neighboring 
particles are allowed to leave or enter the system in ordered groups: 

ABC +± 000, (8) 

p 

where /x is a chemical potential, taken to be equal for all three species and p is a 
parameter whose value does not affect the steady state in the case where detailed balance 
is satisfied. This particular form of the nonconserving process is chosen so that the equal 
densities condition, Na = Nb = Nc = N/3, could be maintained. 

In general, this model, consisting of the dynamical rates ([I]),© and (jSJ), exhibits 
nonvanishing currents in the steady state, similarly to the original ABC model. 
However, as demonstrated below, for equal densities the model exhibits detailed balance, 
reaching an equilibrium state with a distribution: 

P{{X i }) = ±-q H <*>« x *», (9) 

governed by the Hamiltonian: 

U GC ({X}) = H ({X}) - \n (N - 1) - /iNL. (10) 

o 

Here H ({Xj}) is the Hamiltonian of the standard ABC model, as given in Eqs. (j3J) or 

(ED. 
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We now verify that the generalized model indeed obeys detailed balance with 
respect to Hoc- Under particle-conserving processes ([1]) and ([?]), detailed balance is 
maintained due to the fact that vacancies do not affect the energy of a configuration. 
For the nonconserving process (jSJ), detailed balance is verified by noting that the 
energy of a configuration is invariant under translation of ABC triplets. Namely, 
E(. . . Y ABC • • •) = E(. . . ABCY . . .), where E is the energy and Y stands for either 
a particle of any species or a vacancy. Thus, the change in energy due to depositing 
or evaporating ABC triplets is independent of where on the lattice this process takes 
place. Depositing a triplet of ABC into a 000 triplet changes the total energy (|T0|) of a 
system with an initial particle number N by 

AU GC = N + 1 - \ (6A + 6) - 3/iL = -3/iL, (11) 
o 

which is consistent with the local dynamical rates in Eq. ([H]). Therefore, detailed 
balance is maintained for the nonconserving process as well. Equation (TTTj) implies that 
the dynamical parameter fi is in fact the conjugate field of N. 

We proceed by considering two cases. The first is a conserving model whose 
dynamical rules consist of Eqs. (CQ) and (J?]). The second is a nonconserving model 
that evolves by all three types of processes given in Eqs. ([I]), ([7]) and ([H]). With the 
effective Hamiltonian, l-icc-, the two types of dynamics correspond to the canonical and 
the grand canonical descriptions of the ABC model, respectively. 

4. Phase diagram of the ABC model with conserving dynamics 

We consider the generalized ABC model under conserving dynamics, in the case of 
equal densities and in the weak asymmetry limit, q = e~^l L . Previous studies of the 
standard ABC model for equal densities [32] found a second order phase transition from 
a homogeneous to an ordered phase at f3 c = 27ry / 3- This result can be easily extended 
to the generalized model with conserving dynamics by noting that the vacancies do 
not contribute to the energy of the model, and thus they are randomly spread in the 
lattice in the equilibrium state. This allows us to map each microscopic configuration of 
the generalized model to that of the standard ABC model by removing the vacancies, 
yielding a 'condensed' system of size N. The mapping may be reversed by adding 
to the 'condensed' system the L — N vacancies in all possible arrangements with equal 
probability, resulting in a one-to-many correspondence. The free energies of system with 
L — N vacancies and its 'condensed' counterpart, denoted as F(Na, Nb, Nq, L — N) and 
F(Na, Nb, Nc-i 0), respectively, differ only by a shift due to the entropy of the vacancies: 

Here and throughout this paper the free energy is rescaled by j3. The 'condensed' N- 
size system is thus equivalent to the standard ABC model. Writing q = exp (—/3/L) = 
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exp (— 0r/N), where r = N/L, it is evident that the 'condensed' system has an effective 
inverse temperature of 0r. It therefore exhibits a phase transition at the critical line 



0c 



(13) 



which is also the transition line of the corresponding generalized ABC model of length 
L. For r = 1, the model contains no vacancies, and we recover the transition point of 
original ABC model. 

In order to compare the phase diagrams of the conserving and the nonconserving 
dynamics, we plot the phase diagrams in the (1/(3, p) plane. The chemical potential of 
the conserving model is obtained using the Hamiltonian 



'Hc({X i })=H({X i })--N(N-l) 

o 



(14) 



with H ({Xj}) given by Eqs. (jlj) or (jSJ). The last term in Tic is chosen so that the 
Hamiltonian differs from the nonconserving Hamiltonian (TTOl) only by the term —pNL, 
as required when defining the canonical and grand-canonical Hamiltonians of a model. 
The free energy of the conserving model in the homogeneous phase is given in the large 
N and L limit by 



N N N Ar 

■ rA rrr L ~ N 



AH n 



N) In (L — N). 



(15) 



Due to the specific choice of the constant energy shift in the Hamiltonian (JHJ), the 
energy vanishes in the homogeneous state, and only the entropy contributes to the free 
energy. The chemical potential in the homogeneous phase is thus given by 



la? 7 * 



dN / 
and the critical line, where r 



In 



N 

T )-ln(L-A0 



27r\/3//5, can be written as 



In 



2vr 



In 1 



2tt\/3 N 



(16) 



(17) 



The resulting phase diagram of the conserving model is shown in Figure [TJ 

We now provide an alternative derivation of this phase diagram by expanding the 
free energy of the model in small deviations of the density profile from the homogeneous 
solution. This approach yields more information about the nature of the transition, 
and will especially be useful in the next section for analyzing the nonconserving phase 
diagram. 

For this purpose, we turn to the continuum limit 32j, |36|], where the local densities 
of A, B and C particles at the point x = i/L are represented by the density profile 
p n (x), n = A,B,C, with p(x) = Pa( x ) + Pb(%) + Pc( x )- The average particle density 
is r = N/L = J Q p(x)dx. The steady-state distribution of the density profiles may 
be expressed as P [p n (x)\ = exp {— LJ 7 [p n (x)}}, where jF[p n (x)] is the free energy 
functional, rescaled by 0. The equilibrium profile can thus be found by minimizing 
the free energy functional with respect to p n {x), under the equal densities condition, 
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Figure 1. The (/^, phase diagram of the generalized ABC model with conserving 
dynamics, displaying a second order phase transition line between a homogeneous and 
a phase-separated states. Typical density profile, p n {x), in each phase are given in the 
insets. 



Jq pa{x)<Ix = Jq pb(x)cIx = f pc{x)dx = r/3. The free energy functional for the 
conserving model is 

F[p n (xj\ = dx [pa(x) ^Pa(x) + p B (x) lnp^x) 



+ p c (x) \npc(x) + (1 - p(x)) In (1 - p(x))) 



+ dx J dz[pA{x)p B (x + z) + pb(x)pc(x + z) 

+ pc{x)p A {x + z)}z- ^r 2 \, (18) 

where the first integral corresponds to the entropy and the second integral corresponds 
to the continuum limit of the Hamiltonian (TH1) . 

Some characteristics of the equilibrium profile p n (x) can be determined by the 
symmetry of the model. Due to the cyclic boundary conditions, Pa{x), Pb{%) and Pc{%) 
are periodic functions with period 1. In addition, since the Hamiltonian favors phase 
separation between the three species, we expect the density profiles in the equal densities 
case to satisfy 

Pb(x) = Pa (x - , p c {x) = p A (x+ M . (19) 



This assumption is further justified in |Appendix C[ We now use the translation 
symmetry of the model and set x = at the symmetry axis of Pa{x)- The coarse 
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grained state of the model can thus be represented by a Fourier series for Pa(x) 

oo 

- + a m cos (2m7rx) . (20) 



p A {x) 



m=l 



The relation (Tl9|) implies that the other two profiles are given by 

r \ -< / 1" 

Pb{x) = - + 



am COS 



Pc{x) 



m=l 



2rmr [ x — — 



(Xffi COS 



m=l 



2m7T ( x + 



At high temperatures, T = all coefficients a r 



(21) 

vanish, and the profile is 



homogeneous, with pa(x) = pb(x) = pc(x) — r/3. 

In order to find the transition line between the disordered and ordered phases we 
expand J 7 close to the homogeneous profile in terms of a small perturbation by assuming 
a m <C 1. The amplitudes evolve by ^f- = — In |Appendix A| we show that the 



dt da, 

transition to the inhomogeneous phase takes place at (3 



2n\/3/r when the first mode, 
Oi, becomes unstable, whereas all higher order modes are linearly stable. Just below this 
critical line the higher order modes (m > 1) are driven by a\ and may be represented 
by a power series of a\. The amplitude of a\ may thus serve as the order parameter of 
the transition. 

The amplitudes of the higher order modes (m > 1) are obtained by setting 
dF/dcim = 0, which yields to lowest order a m ~ a™ (see Appendix A ). The fact 
that the vacancies are homogeneously distributed in the equilibrium state implies that 
the total particle density is constant in space, po(x) = 1 — p(x) = 1 — r, and hence 
5pa{x) + Spb{x) + Spc{x) = 0. From this it follows that all a^, m coefficients vanish for 
m > 1. Consequently, the expansion of J-"[p ra (a;)] in powers of a\ is greatly simplified. 
Up to order a\ it requires terms that involve only the coefficients a x and a 2 , yielding 

3V/3/A 2 
8tt J 1 



F[r/3 + 6p n (x)} = J^ h (r 



+ 



9 3^3/3 ^ 




a 2 - —0^2 + 



if + O (a*), 



Ar ' 16vr I " z 8r 2 ~ L ~~* ' 32r 3 
where J r h( r ) is the free energy of the homogeneous profile given by 
J^h ( r ) = t In + (1 — r) In (1 — r) . 



(22) 



(23) 



Details of this derivation are given in Appendix A[ 

We can express a 2 in terms of ai using the equation dJ- ' jda^ = 0, which yields: 

9tt 



a 2 



(24) 



r (VSf3r + 12tt) X * 

We finally obtain the following Landau expansion of the model, given by the power 
series of J 7 in the order parameter a\\ 

7 [p n {x)\ = F h {r) + f 2 a\ + Ua\ + O (a\) , (25) 
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where 



f2 (ft r) 



9 

4r 



3\/3/0 



57T 



/ 4 (/3,r) 



81 / V3[3r + 6n 



32r 3 



12tt 



> 0. 



(26) 



y/30r 

By setting f 2 = we obtain the same critical line as in Eq. ( TT3l) . shown in Figure [U The 
fact that / 4 is positive for any value of r indicates that this is a second order transition 
line, from the disordered phase, where f 2 > and T is minimized by a\ = 0, to the 
ordered phase, where f 2 < 0, in which J 7 is minimized by a nonvanishing ai. 



5. Phase diagram of the ABC model with nonconserving dynamics 

5.1. The second order line 

The free energy functional corresponding to the generalized ABC model with 
nonconserving dynamics is 

g[pn{x)}=J r [p n {x)}-(3fir, (27) 

where jF[p n (x)] is given by Eq. ( TTBl) . In order to find the transition between the 
disordered and the ordered phases, Q may be expanded close to the homogeneous profile, 
as was done in the previous section in analyzing J 7 . Here, however, in addition to the 
modulation of the density profiles of A, B and C, parameterized by a m , deviations of the 
overall density, r, represented by Sr, must be taken into account. Thus, the A-particle 
density profile close to the transition can be written as 



Pa{x) 



r- OO 

r or sr^ 



cos(27rmx) 



(28) 



m=l 



where here, again, ps(x) = pa{x — 1/3) and pcipc) = Pa(% + 1/3). Similarly to the 
conserving model, the flat profile of the vacancies implies that all a^ m = for m > 1. 
One can show that to lowest order a m ~ a™ and 5r ~ a\. This simplifies the expansion 
of Q considerably. We carried out the expansion to eighth order in a±, however in order 



to avoid lengthy expressions, we outline it here and in Appendix B only to sixth order 



In Appendix B we find that 



9 

81 



+ 



2(1 -r) 2r 



9 3^3/^ 



4r ^ 16ir 



{Sr 



, 2 27 



yr 



a?a 2 



2"1 



, . 4 243 6 9 2r 81 2 2 

+ 32r^ ai + 32^ ~ 4^ A + S^ 2 
9 2 2 243 4 243 4 

4r 3 1 32r 4 1 16r 4 1 

27 2 r f 1 

+ — ala 2 dr + — 

4r 3 V6 1-r 2 



(29) 
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g h (r) = r In ( T - ) + (1 - r) In (1 - r) - /3/ir (30) 



where 

(VI = r In I 

,3 

is the free energy of the homogeneous profile. 

In order to characterize the nature of the transition line of the nonconserving model, 
the expansion in Eq. (l29p has to be continued to eighth order in ai, by taking into 
account terms that involve only the amplitudes and Sr. The amplitudes are 

substituted by the following series in a\. 

5r = A 0)2 al + A 0A a\ + A Qfi a\ 
a 2 = A 2 ^a\ + A 2A a\ + A 2fi a\ 

04 = A^a-y. (31) 
The coefficients {Ajj} are derived from the equilibrium condition dQ/d(5r) = and 



gi ((3 c ,r) = —(3r-l) . (34) 



dQ/da m = for m > 1 (see Appendix B). Thus, the Landau expansion of Q in terms 
of a\ is obtained: 

Q [p n {x)\ = Q{r)+ g 2 a\ + g A a\ + g & a\ + g 8 af . . . (32) 
The second order coefficient 

fc(0,r) = / a (Ar) = f (33) 

4r 8tt 

vanishes at (3 C = 2n\^3/r. On the critical line, f3 = (5 C) the fourth order coefficient is: 

27 
32~7 3 

It is positive for r > 1/3 and it becomes negative for r < 1/3. Therefore, there is a 
multicritical point (MCP) at tmcp = 1/3, with 

/W = Sttv^ ~ 32.648 /xmcp = ~ -0.0549. (35) 

67rv3 

Calculating the higher order coefficients, we find that on the critical line 

S 6 (/3 c ,r) = ^(6r 2 -5r + l), (36) 
27 

98 (/?c ' r) = 1024^ ( 1215r3 ~ 1692r2 + 762r ~ 109 ^ ' (37) 
The sixth order coefficient vanishes at the MCP, g§ (/9 c ,^mcp) = 0. This seems to 
be an accidental coincidence, and it will be discussed in more detail in Section [6j The 
coefficient becomes negative just above the MCP. However, the eighth order coefficient 
is positive, g$ (Pc^mcp) > 0, and large enough so that the second order transition is 
stable above the MCP. Thus, the MCP is, in fact, a fourth order critical point. A plot of 
the coefficients g±, g§ and #8 along the critical line as a function of r is shown in Figure 
[2j At low densities, r < Tmcp — 1/3, #4 (/3 C , r) is negative and we expect the transition 
to become first order. 
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Figure 2. The 4 th fi th and 8 th order coefficients in the expansion of Q along the critical 
line. Both the 4 th and 8 th order coefficients are positive above the MCP. The 6 th order 
coefficient is negative just above the MCP, and vanishes at the MCP itself. 



5.2. The first order line 

In this section we complete the phase diagram of the nonconserving model by evaluating 
the first order transition below the multi-critical point. We begin by considering the 
behavior of the model at 1//3 — T — 0, where the density profile has a trivial form. 
We then derive the first order line using an analytic expression of the mean-field profile. 
Finally we provide a simple approximation for the first order line at low temperatures, 
which also yields a lower bound for the transition for arbitrary temperature. 

In the limit = T — > 0, the entropy may be neglected and the free energy 
functional, Q[p n ( x )], is given by the energy: 

Gt=o \Pn( x )} = P jjf dx j dz[p A (x)p B (x + z) + p B (x)pc(x + z) 

+Pc(x)p A (x + z)]z -^r 2 - //rj . (38) 
It is straightforward to verify that the ground state profile is the fully separated state, 

pY(x) = I r f < 5 , (39) 
HA V ! 10 otherwise ' V ; 



with p s B P (x) = p s A P (x — 1/3) and p^ p (a;) = p a ^(x + 1/3). The free energy of this state, 

18 



g T = (r) = -f3[^-r 2 + pr), (40) 
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is minimal at r = for p < — and at r = 1 for /i > — -L. Consequently, there is a 
discontinuous transition from an empty system to a fully occupied phase-separated one 
at fi = — and T = 0. This suggests that there is a first order transition line, denoted 
here as (T*,p*), which connects (0, — i) and the MCP. 

The derivation of the first order transition line requires an explicit expression for 
the density profile, p n {x), at finite temperatures. To this end we first compute the 
profile of the conserving model using its mapping to the standard ABC model presented 
in the first paragraph of Section |H According to this mapping, the steady state of 
a conserving model of size L can be extracted from that of the standard ABC model 
(without vacancies) of size N with an effective inverse temperature of fir. For the 
standard ABC model we can apply the analytic solution of the mean-field equations 
which has been suggested by Fayolle et al. and derived explicitly by Ayyer et 



al. [36] (see also Appendix C ). The solution has been formulated for the ABC model 
on an interval, but for the case of equal densities it applies also for periodic boundary 
conditions. We use it to obtain the profile at an inverse temperature /3r, and map it back 
to the profile of the corresponding generalized model (with vacancies) by multiplying it 
by r, yielding : 

. 1 + sn (2f3rx/x, k) 

PA(x) = r ' \ , (41) 

a+ — a_sn {zfirx/ k) 

where sn stands for the Jacobi elliptic function, and x, a + ,a_,k are functions of the 

parameter fir whose form is given in Appendix C The profiles for B and C are again 

given as Pb(x) = p A (x - 1/3) and pc(x) = p A (x + 1/3). 

As shown in Appendix C the density profile in Eq. (T4TT) is also a stationary solution 

of the mean-field equations of the nonconserving model. These mean-field equations 

include, however, an additional constraint, 

pi ( x ) = e ~ 3 ^p A (x) p B (x) p c (x) , (42) 

which results from the detailed balance condition relating the evaporation and deposition 
processes OH]). This constraint yields the relation between p,r and given by 



rK 1 ' 3 {fir) 



(43) 



1 — r 

where K is independent of x, and obeys Kr 3 = p A {x) ps{x) pc{x) . The dependence of 
K on the parameter fir is given in Appendix C Equation (14*31) also defines the chemical 
potential in the conserving model, where each steady-state profile is also a stationary 
solution of the nonconserving model with that value of p. 

The resulting p (r) curves for fixed fi are shown in Figure [3j The key feature is the 
region of p in Figure[3b where there are 3 available solutions for r. The solution with the 
intermediate values of r has negative compressibility and it is therefore unstable under 
the nonconserving dynamics. At the value of p for which the two other stable solutions 
have the same free energy (denoted by a dashed line) the nonconserving model undergoes 
a first order phase transition. The transition point, p*, is found by evaluating the free 
energy of the nonconserving model, Q, using the chemical potential and the density 
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Figure 3. The /i(r) curve at constant temperature as calculated in the conserving 
model using the analytic solution of the mean-field equations. The thin and thick lines 
correspond to the homogeneous and ordered phases, respectively, and the transition 
point is marked by x. In (a) T = 0.04 > T MC p, while in (b) T = 0.02 < T M CP- The 
dashed line in (b) marks the point where the free energy of the two phases is equal 
and the nonconserving model exhibits a first order transition. 



profile given above. In principle, one can deduce the entire phase diagram using this 
nonperturbative approach. However, the expansion of the free energy in Section 15.11 is 
more convenient for characterizing the nature of the transition. 

The phase diagram of the nonconserving model with the resulting first order line is 
presented in Figure HI The numerical evaluation of Eq. ( l4Tj) and ( l43|) requires numerical 
precision that grows linearly with — 1/T making it prohibitive at low temperatures. 
To study this limit we now present a simple approximation for the transition line which 
also yields a lower bound for the transition at all temperatures. 

We first observe that above the transition curve (T > T* or fi < /i*) the optimal 
profile is homogeneous, with p n (x) = r/3. The value of the average density r is 
determined by minimizing the free energy of the homogeneous profile (130|) with respect 
to r. As the model approaches the transition point, the free energy, Q, develops a local 
minimum corresponding to a inhomogeneous density profile. The transition temperature 
may be defined as the lowest temperature for which the free energy of the homogeneous 
profile is lower than that of any other profile. Instead of considering any possible profile, 
we look at the totally separated profile in Eq. (139|) . whose corresponding free energy is 
given by 



For T > the model eventually displays an ordered phase that is not fully separated 
and therefore has a lower free energy (due to entropic effects) than Q scp (r). Hence, for 
a given \x the temperature at which 




(44) 




(45) 
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Figure 4. Phase diagram for the generalized ABC model with nonconserving 
dynamics. The second order (thin line) and first order (thick line) transitions are 
separated by a multicritical point (★). The first order line terminates at (T — 0, /i = 
— 1/18). The inset displays the lower bound for the first order line (dashed line) in 
comparison with the exact first order transition (see text). 

is a lower bound for the transition temperature, T*. At low temperatures, the first order 
transition line approaches this lower bound and the ordered profile is well represented 
by the fully phase-separated one. Thus, this bound provides a good approximation for 
the behavior of the system at low temperatures (T <C Tmcp), as shown in the inset of 
Figure HI 

5.3. Monte Carlo simulations 

The picture emerging from the continuum limit is supported by the results of Monte 
Carlo (MC) simulations, performed under effective equilibrium conditions, = Nb = 
Nq = N/3, with both conserving and nonconserving dynamics. 

In the nonconserving simulations, the state of the system is updated according 
to the following procedure: At each time-step a single site is randomly selected. The 
type of move to attempt is also chosen randomly: with probability 1/2 an attempt is 
made to exchange the particle of the chosen site with its right-hand neighbor and with 
probability 1/2 an evaporation or deposition process is attempted. In the latter case, 
if the chosen site is occupied by a B particle with an A to its left and a C to its right, 
then the removal of the triplet ABC is attempted. If the chosen site is vacant, and it is 
surrounded on both sides by vacant sites as well, the condensation of an ABC triplet is 
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attempted. Finally, the move is accepted with probability 

exp(-d^) A/-/.,,- I. ^ 
1 otherwise ' 

where AHgc is the change in energy due to the chosen move, as determined by Eq. ffTUl) 
(possible values are 0, ±1 and ±3fiL). This procedure leads to a steady state which 
obeys detailed balance with respect to the nonconserving Hamiltonian (TTU]) . 

In order to compare the results of these simulations with those of the conserving 
simulations one has to calculate the chemical potential, /i, for a given temperature under 
conserving dynamics. This can be done by employing in the conserving simulations a 



method similar to the Creutz algorithm for microcanonical MC simulations [39]. The 
idea is to perform a constrained nonconserving simulation so that the average density, 
r, is allowed to fluctuate only below its initial value, r(0), while extracting the value 
of \i from these fluctuations. The simulation is executed as in the nonconserving case 
by selecting a site and an attempted move. Steps in which neighboring particles are 
exchanged are accepted according to Eq. (H6|) . The particle nonconserving steps are 
performed in conjunction with an additional single degree of freedom, termed 'demon', 
that exchanges particles with the system. The 'demon' is initially empty. An attempt 
of removing ABC triplet is accepted with probability 1, and the removed particles are 
added to the 'demon'. Steps that require the deposition of an ABC triplet on three 
vacant sites are accepted only if the 'demon' is not vacant and with a rate given by Eq. 
( 1461) with 1-Lgc replaced by the canonical Hamiltonian, He ( H^l) . The triplet is then 
removed from the 'demon'. As a result the average density is allowed to fluctuate, but 
only to states with an average density below r(0). Fluctuations with higher densities are 
rejected. In the thermodynamic limit this procedure yields the canonical distribution 
of the system with density r(0), even for systems with negative compressibility [20]. 
The probability distribution of the number of particles in the 'demon', P(Ndemon), is 
recorded during the simulation. The chemical potential, fi, is determined from this 
distribution using the equilibrium relation P (iV demon) ~ exp (— f3fiNd em0 n)- A typical 
distribution is given in Figure [5[ from which /i is extracted by applying a linear fit to 
ln[P (N demon )}. 

The simulation results are displayed and compared with the mean-field solution 
(|43|) in Figure (6[ where the average density, r, is plotted as a function of the chemical 
potential, fi. Above the multicritical point, in Figure Ek, we see that the conserving and 
nonconserving simulations follow the same curve. The two types of dynamics are thus 
equivalent, and undergo a second order transition at the critical point (fT7l) marked in 
the figure by x. Below the multicritical point, in Figure [6b, the conserving simulation 
exhibits a similar second order transition, whereas the nonconserving simulation shows a 
discontinuity in r. At the intermediate-density states we find negative compressibility in 
the conserving simulation. These states are unstable under the nonconserving dynamics, 
where r can fluctuate freely. The nonconserving model thus exhibits a discontinuity in 
r, accompanied by hysteretic behavior. This hysteresis is an indication of the first order 
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Figure 5. Natural logarithm of the probability distribution of the 'demon' occupation 
number as obtained by MC simulations of the conserving dynamics for L = 1800, 
N = 1125 and j3 = 25. Linear fit (dashed line) of the simulation points yields /i ~ 0.047. 
Such simulations yield the fi(r) curve. 



transition. The value of /i at the transition, as found in the thermodynamic limit using 
the analytical procedure discussed in Section 15.21 is denoted by the dashed line. This 
thus provides a direct observation of the inequivalence of two ensembles. The results of 
the simulations fit very well the mean-field solution (143|) . 

6. Further generalization of the model 

The nonconserving ABC model with equal densities exhibits a seemingly accidental 
coincidence whereby the sixth order coefficient in the expansion of the free energy 
vanishes at the MCP, i.e. gi = g± = g$ = 0. If indeed this feature is accidental 
one expects any modification of the model to remove this degeneracy. In this section we 
consider a simple generalization of the nonconserving dynamics that maintains detailed 
balance under the equal densities condition. Thus, the features of the phase diagram 
may be found using the expansion of free energy functionals, as presented in the previous 
sections. 

The generalization consists of replacing the factor 1/6 in the Hamiltonian in Eq. 
f fTUj) by a free parameter 7, 

U GCa ({Xi}) = U ({X t }) - 7 iV (N - 1) - fiNL, (47) 

where % ({Xi}) is given by Eqs. ([41 or ([3]) . By imposing the condition for detailed 
balance with respect to the distribution P 7 ({Xj}) = q^ooMi^i) JZ^ for systems with 
equal densities Na = Nb = Nq, and reversing the argument of Eq. (1111) . we find that 
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Figure 6. The fj,(r) curve of the conserving (A) and nonconserving (•) dynamics 
obtained by direct simulations for L = 1800. The simulations were performed at two 
temperatures (a) T = 0.04, above the MCP, and (b) T = 0.02, below the MCP (T MC p ~ 
0.03). In (a) both the conserving (A) and nonconserving (•) simulations result in the 
same curve within the numerical accuracy. In (b) the two types of dynamics yield 
different curves, with the nonconserving dynamics displaying a discontinuity in r at a 
first order transition. Hysteretic behavior is observed in the nonconserving case. The 
expected first order transition point, obtained by minimization of Q, is indicated by 
the dashed vertical line. The solid lines correspond to the analytic solution of the 
mean field equations and so is the x, which denotes the second order transition point 
in conserving model. 



detailed balance is maintained for the following evaporation and deposition rates: 



pq <-»<-,, 7 



ABC +± 000, (48) 
v 

where 

An G c, 7 = 3/iL - (N + 1) (1 - 67) • (49) 
Thus the non conserving dynamics corresponding to the modified Hamiltonian (T4j 



consists of the processes (JTJ), (J2J) and ( I48jl . For 7 ^ 1/6 the evaporation rate depends 
on the particle number, N, and therefore the dynamics is nonlocal. 

We now analyze the phase diagram corresponding to this generalized model in the 
limit of weak asymmetry, q = exp (—j3/L). The continuum- limit calculation of the phase 
diagram can be repeated, producing Landau expansions of the free energy functionals J-" 7 
and Qj. In the conserving dynamics the particle number iV is a constant, and hence the 
7-term has no effect on the expansion of J-" 7 . Therefore, regardless of 7, the conserving 
model exhibits a second order transition at /3 C = 2ir\/3/r. 

The expansion of Q 1 , as detailed in Appendix B[ is given by, 



£7 \Pn{x)] = £ 7 (r) + g\a\ + g\a\ + gja* + g^af + (50) 
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where g% = #2 is independent of 7. It vanishes on the critical line (3 = (3 C = 2iT\/3/r. 
On this line the fourth order coefficient is given by 

9r - 3 + 4nV3 (1 - 67) (1 - r) 
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3 + 27iV3(l- 67) (1 



(51) 



which is positive at all densities satisfying 

V3-Ml-6 7 ) 

' > ' M0P = 3^5 -4„ (1-67)- m 
The coefficient vanishes at r = tmcp which yields: 

6tt [3V3-47r(l-67)] 
PMCP ~ 3 - 4tt V3(l - 67) 

A^mcp = \ [3 - 47iV3(l - 67) 
o L 

2(1-67) , ln^(3-47ry^(l-67))] ] 
9 - 4ttV3(1 - 67) 7T (3\/3 - 4tt(1 - 67)) J 

Note that as 7 — y (47r — v3) /247T ~ 0.144 (from above), the density tmcp vanishes, 
and #4 is positive for all values of < r < 1. 

The phase diagram in the vicinity of the MCP depends strongly on the sign of the 
sixth order coefficient, gj, at the MCP where gj = gj = 0. From the expansion of 



(53) 



we find 

, 7 VSn (1 - 67) [3 - 2V37T (1 - 67)] [4737T (1 - 67) - 9] 5 

(ife) MCP — r F ', TT5 ' ( 5 > 

8 [3 - 4^ (1 - 67)] 

For 7 > |, one has (<7g) M cp > 0> an d the MCP is a tricritical point (TCP). Thus, 
the phase diagram in this case consists of a second order line which becomes first order 
below the TCP. However for 7 < |, the coefficient (g2) MCP is negative and the trictitical 
point is unstable. As a result the first order line intersects the second order line above 
the TCP. Thus the phase diagram consists of a critical line which terminates at a first 
order line at a critical end point (CEP) where g\ = 0, gj > and g$ < 0. This point 
is located above the TCP (which is unstable). The first order line continues into the 
ordered phase, and it ends at a critical point (CP), as shown schematically in Figure [71 
Above the CEP the first order line marks a transition between two ordered phases with 
low (I) and high (II) densities. As 7 approaches 1/6 from below, the CEP, CP and TCP 
approach each other. They merge at 7 = 1/6, yielding a fourth order critical point. 

To complete the analysis of the phase diagram one has to determine the first order 
line. This can be done by using the exact solution for the density profiles, as discussed 
in Section 15.21 The location of the first order transition at T — y can be derived by 
equating the free energies of the homogeneous and the fully phase-separated states 

min£ 7ih (r) = min £ 7jSep (r), (55) 
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Figure 7. Schematic phase diagrams of the modified nonconserving ABC model with 
nonlocal dynamics for (a) 7 > h and (b) 7 < ^. In both cases the phase diagram is 
composed of a second order transition (thin line) at high temperatures, which becomes 
first order (thick line) at low temperatures. The dashed line marks the continuation of 
the second order line, which is unstable for the nonconserving dynamics. The first order 
line intersects the T = axis at fi = ^ — 7. In (a) The two lines meet at a tricritical 
point (TCP), as in the local dynamics case, 7 = i. On the other hand in (b) the 
first order transition line intersects the second order line at a critical end point (CEP). 
The second order line terminates at this point while the first order line continues into 
the ordered phase, where it marks a transition between two distinct phase-separated 
states, with low (I) and high (II) densities, and terminates at a critical point (CP). The 
hollow point on the unstable second order line, denotes the TCP where 52 = .94 = 
and gg < 0. This point and the dashed segment of the second order line are preempted 
by the first order line and are thus not accessible within the nonconserving dynamics. 
For particles conserving dynamics the phase diagram exhibits a second order line, 
composed of the thin and dashed lines in the figure, for any value of 7. 



where 

7 ,h(r)=rln(0 + (l 

and 

£ 7 ,sc P (0 = rln (r) + (1 

At T = the model exhibits a discontinuity in the total density from r = to r = 1 at 
fi = 1/9 — 7, and thus undergoes a first order transition at that point. 

Figure [7] shows a schematic representation of the phase diagrams of the generalized 
ABC model with nonlocal dynamics for conserving and nonconserving dynamics above 
and below 7=1/6. As expected, the nongeneric feature of a fourth order critical 
point has been removed by slight modification of the model. Generically, we expect the 
multicritical point to become either a tricritical point or a critical end point, as shown 
in the figure. 



r) In (1 — r) + /3 ( — — 7 ) r 2 — /3/ir, 



(56) 



r) In (1 — r) + /3 ( — — 7 ) r 2 — (3 fir. 



(57) 
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7. Conclusions 

In this paper we generalized the ABC model to include vacancies and processes that 
do not conserve the particle number. This enables us to analyze and compare the 
phase diagrams of the conserving and the nonconserving models. We have shown that 
in the case where the average densities of the three species are equal, the dynamics 
of the generalized model obeys detailed balance with respect to a Hamiltonian with 
long-range interactions, despite the fact that the dynamics is local. Studying the (/i, T) 
phase diagrams of the model for equal densities, we found that in the conserving case 
it is composed of a second order line separating the homogeneous and phase-separated 
states, while in the nonconserving case the second order line becomes first order at low 
temperatures. The analysis of the phase diagram has been carried out by studying 
the Helmholtz and Gibbs free energies for the conserving and nonconserving dynamics, 
respectively, in the continuum limit. As has been shown in the past this limit yields 
the exact steady states in the thermodynamic limit. In this study we applied a critical 
expansion of the free energy near the homogeneous phase and an exact solution of the 
mean-field equations for the density profiles of the phase-separated state. The results of 
this analysis are verified by direct Monte Carlo simulations of the two types of dynamics. 

The fact that the two types of dynamics result in rather different phase diagrams 
can be associated with the inequivalence of the canonical (conserving) and grand- 
canonical (nonconserving) ensembles in systems with long-range interactions. We find 
that as expected from studies of long-range interacting systems, the two ensembles yield 
different steady states in the region where the grand canonical ensemble displays a first 
order transition. 

We expect the generalized ABC model to display similar behavior even for small 
deviations from the equal densities case, where detailed balance is not satisfied. The 
present study can thus serve as a starting point for a study of the ABC model out of 
equilibrium. It could provide an interesting correspondence between some properties of 
the well-understood equilibrium systems with long-range interactions and those of the 
less well- understood nonequilibrium driven models. Details of a study of the generalized 



ABC model with unequal densities will be published elsewhere 40|. An interesting 



driven model where ensemble-inequivalence has recently been observed is the zero- 



range process 
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In this model, the drive, provided by the spatial asymmetry of 
the transition rates, does not influence the steady state. This state thus remains the 
same as the steady state of the nondriven equilibrium model, where the transition rates 
are symmetric, which can be expressed in terms of a Hamiltonian. By contrast, the 
ABC model provides a framework within which one can readily probe nonequilibrium 
steady states which are not expressed by a Hamiltonian. 
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Appendix A. Critical expansion of the conserving free energy 



We present the critical expansion of the free energy of the conserving model. In the 
continumm limit, the free energy, rescaled by /3, is given by T = fie — s, where 



s[p n (x)] = - dx [p A (x)\n(p A (x)) + p B (x)]n(p B (x)) 
Jo 

+pc(x) In (pc(x)) + (1 - p(x)) In (1 - p(x))} (A.l) 

is the entropy per site of the profile, derived from simple combinatorial considerations, 
and 



1 f 1 r 

dx dz[pA{x)p B {x + z) + Pb(x)pc(x + z) 



o Jo 



1 



+Pc(x)pa(x +z)]z- ^r 2 

o 



(A.2) 



is the energy per site given continuum limit of the conserving Hamiltonian (IM|) . for 
n = A,B or C. 

We begin by expanding the steady-state profile close the homogeneous solution in 
the most general form : 



Pn(x) = ~ + Sp n (x) = - + a n, m e 



2mnix 



(A.3) 



m=— oo 



For the profile to be real we assume a n - m = u* nm . The time evolution of these modes 



is set by 



(la n 



dt 



da 



— — which yields to lowest order in a n , 



d_ 

If 



o, 



i/3 
2irm 



( 



V 



i 
o 



i 





\ 



/ 



a. 



(A.4) 



where a m = {atA,mi a B,m-i a C,m) an d I 1S the 3x3 unit matrix. For the m mode, the 
highest eigenvalue of the matrix above is p^Jt _ 3 w j^ n ^ ne eigenvector (1, e -2 "" 1 / 3 , e 2?rl//3 ). 
As the temperature is decreased (/3 increased) the first mode to become unstable is 
m = 1 (at /3 = 2n\/3/r), while the other modes are linearly stable. Just below this 
transition line the higher modes (m > 1) are driven by the m = 1 mode. The m th mode 
is driven to lowest order by a term of the form a™ x a n - m in the Taylor expansion of the 
logarithmic function in the entropy. We thus obtain to lowest order a n ^ m ~ a™i with 
the eigenvector (1, e - 2mm / 3 ^ e 2mm/3y ^y e can therefore simplify our expansion and set 



Pa{x) 



+ 5pa{x) = - + am cos(2mirx) 



(A.5) 



771=1 
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with Pb(x) = Pa(x — 1/3) and pc(x) = Pa(x + 1/3)- This latter form is also justified in 
the work of Ayyer et al. |36|, who proved that the ordered profile of the model is unique 
and obeys this symmetry. 

As explained in Section HI the flat profile of the vacancies implies that 5po = 
— ($Pa + 5pB + 8pc) = and therefore a 3m = 0. We wish to evaluate the free energy up 
to order a\ and thus expand the entropy and energy only in terms of a± and a 2 . 

The entropy of the A particles is given by the first term in the RHS of Eq. (|A.1|) : 



sa 



dx 



dx 



+ <5p A (x) 



In 



L3 



+ 6pa(x) 



r (r 
- In - 
3 V3 



+ 



1 + In 



3/J 



3 9 

— (5p A (x)f + — {8p A {x)t 



After performing the integral we obtain: 



r. fr\ 2 
— Sa — — In — + Oi 
A 3 \3/ 1 



4r 



4 ( 27 \ 2 
1 I 32r 3 J 1 



+ «2 



9 



3 2 
-0(5p A (x)f. 



3 

4r 



(A.6) 



(A.7) 



Under the condition of equal densities one has sa = Sb = sc, and therefore the total 
entropy of the particles is equal to 3sa- The total particle density is constant in space, 
p(x) = r, and thus the entropy of the vacancies, (1 — r) In (1 — r), contributes a constant 
term to the expansion. The total entropy of the system is thus: 

T\ , / 9 



r In 



+ a 2 2 



9 



1 — r) In (1 — r) + a l 
SI 



a\d2 



4r 
27 



(A. 



4r J " x \32r 3 j ' \tsr* 
Similarly in the expansion of the energy, e, we consider the interaction energy of A 
and B particles given by the first term in the RHS of Eq. (]A.2|) : 



r 2 2 2 

AB 18 8tt 1 16tt 2 
From symmetry, 6ab = £bc — £ca- Adding the constant term — r 2 /6, we find: 

3^ 9 / 3^' 



S7T 



+ a 2 



16tt 



(A.9) 



(A.10) 



The critical expansion of the free energy in Eq. ( 1221) is then obtained by inserting Eqs. 
(IA~8|) and (|A~T0|) into 7 = (3e - s. 



Appendix B. Critical expansion of the nonconserving free energy 



We present the critical expansion of the free energy of the nonconserving model whose 
Hamiltonian (1471) is given by 

V- G c a ({X,}) = H ({X,}) - 1 N (N - 1) - pNL. (B.l) 
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This is a modified Hamiltonian for the case of nonlocal dynamics (see Secition E]). It 
reduces to the Hamiltonian (fTUj) considered in Section [5] by setting 7 = |. The free 
energy of the model is thus 



£7 [P«(»] = 



<-[/>„(•/■)] ( ^ 



7 r 



///• 



s [Pn(a:)] 



(B.2) 



where entropy and energy, s and e, are given in Eqs. flA.ll) and (1A.2|) respectively. 
We follow the same expansion procedure presented in Appendix A where the free 
energy is written as a power series of a\. The analysis of the multi-critical point of 
the nonconserving phase diagram requires terms up to order af. We carried out the 
calculation to this order. However, in order to avoid lengthy expressions, we present 
it here only up to order af, where the free energy is expanded in terms that involve 
only the amplitudes at and a 2 . Here we also take into account fluctuations in the total 
particle density, r, denoted by Sr. The density profile is thus expressed by the Fourier 
expansion, 



p A (x) 



r 5r 

— I h ai cos (2nx) + a 2 cos (Attx) . 

3 3 



(B.3) 



The calculation is similar to that detailed in |Appendix A[ with some modifications. 
First, the dependence of the entropy of the vacancies on the density, and thus on Sr, 
has to be taken in to account. It takes the form: 



so 



+ 



(Sr 



r) In (1 



r) -Sr [1 + ln(l 

\3 



r)] 



(Sr) 



2(l-r) ' 6(1 -rf (B ' 4) 
In addition, the density in the homogeneous steady state depends on the value of \x and 
is determined by the equilibrium condition: 



d_ 

Or 



Gy,h (r) 



In 



3(1 



+ 2/3 ( - - 7 ) r - Pn 



0. 



(B.5) 



where Q 1) \ i is the free energy of the homogeneous profile ( 1561) . The expansion of the free 
energy is: 

' 9 3V3fi 
4r 



£7 \Pn( X )] = ^7,h (r) + 



9 



57T 



4r 



+ 



3\/3/3\ 



9 2, 

+ 



f-^+271 



1 



27 



+ 



rOn + 



243 



32r 3 ^ 1 32r 5 1 
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16r 



4 a x a2 



1 



1 

6r 2 



.6(1 -r) 

We use the following expansion for the amplitudes: 

Sr = A Q)2 a\ + A 0A a* } a 2 = A 2 ^ 2 a\ + A 2A a\. 



(Sr) 3 + O (a?) 



(B.7) 
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Substituting these terms in Q results in the power series: 

£7 KO)] = G 7 ,h 0) + 92 a i + 9~l a \ + gJ a i + (°i) 



where 



and 
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92 



91 
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(B.9) 
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(B.10) 



The coefficient can be expressed in a similar fashion. 

The coefficients are derived from the equilibrium condition: 



dQ 1 \p n {x)\ 







dQ y [p n {x)\ 



0. 



da 2 d (8r) 

Expanding the equation for a 2 in powers of a\ and using Eq. (IB. T|) . we find 

'3^3/3 



(B.ll) 







, -M,2 -5-3 + — A 2i2 ^ 

^7r 8r^ 2r 1 



(-8l7r + 108A 2 , 2 TTr + 36A 0j2 nr - 24A 2:2 A 0i2 irr 2 



167rr 4 

+24A 2i4 vrr 3 + 2^3/3 A 2 , 4 r 4 ) a\ + ... (B.12) 

Each power of a% has to be equal to zero independently. From the second order term, 
we find 

9tt 



.4 



2:2 



r (12tt + \/3/3r) ' 

Similarly, expanding the equation dQ 1 /d{5r) = in powers of ai one finds 

27(1 -r) 



A 



0,2 



(B.13) 



(B.14) 



4r [3 + 0r (1 - r) (1 - 67)] ' 

Equations (IB~13l and ( jBjgJ) are then used to evaluate A 2 4 from the fourth order term 
in flB~T2|) : 

817T 



.4 



2,-l 



X 



2r 3 (I2vr + >/3/3r)' 
v^/Sr 2 [3 + 0(1 - r)(l - 67)] - 18tt(1 - r) 



3 + /3r(l-r)(l-67) 
The higher order coefficients are found in a similar manner. 

Substituting the coefficients in the expression for the free energy Q J} we obtain 



(B.15) 
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3(1 -r) 



V3(3r + 6tt 
V3f3r + 12tt ~ 3 + Jp (1 - r) (1 - 67) 
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To avoid lengthy expressions we display the expression for g 6 only along the critical line 
f3 c = 2ir\/3/r. With the notation 9 = (1 — r) (1 — 67) it is given by: 

kr\ 243 

% (/3c) = — —3 x 

64r 5 [3 + 2tcV36] 

lOvr 3 ^ 3 + 6n 2 9 2 (17r - 5) 

+6nV39 (9r 2 - r - 2) + 9 (6r 2 - 5r + 1)1 . (B.17) 



For 7 = I we find that at the multicritical point, where g\ = gj = 0, one also has 
gl — 0. We therefore need to calculate the eighth-order coefficient, gj, in the same 
manner described above. This requires the evaluation of the amplitude a 4 as well. This 
calculation, whose details are not presented here, yields the following expression for gj 
along the critical line 

kr\ 243 

9i (AO = — ^Ts x 

1024r 7 [3 + 2ny/39] 

5632tt 5 v / 3# 5 + 24vr 4 fl 4 (2883r - 1123) 

+48tt 3 v / 3# 3 (I800r 2 - 717r - 203) 
+72n 2 9 2 (I458r 3 + 567r 2 - 1413r + 268) 
+ 18nV39 (3645r 3 - 3186r 2 + 204r + 217) 

+27 (I215r 3 - 1692r 2 + 762r - 109)] . (B.18) 



The coefficient gj, gj and gj are used for the expansion of the free energy for 7 
Eq. d32]) and for 7 ^ | in Eq. (|50|). 



7; m 



Appendix C. Steady-state profiles in the continuum limit 

In order to locate the first order transition line of the ABC model one has to calculate the 
density profiles of the three species in the ordered phase. In this Appendix we provide 

This 



36 



an analytic solution for the profiles by applying the approach introduced in 
is done first by translating the dynamical rules in Eq. (CQ) to an equation for the time 
evolution of (A4) as 

- (Ai) = q (A-i^) + q (CiA i+1 ) + (B t A l+1 ) + (Ai-iQ) 
at 

-q (AB i+1 ) - q (C^A t ) - (B^Ai) - . (C.l) 

The corresponding equation for (Bi) , (Cj) are obtained by cyclic permutation of A, B 



and C. As has been shown in (36|, in the weak asymmetry limit and for large L one has 

1 



(X l Z l±1 ) = (X l )(Z l±1 ) + 0(- 



(C.2) 



where X and Z are either A, B or C. In the continuum limit 32j, |36|, |37j one can write 
(A i±1 ) =pa± ~t~k~ + od=) (C.3) 
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and similarly for B and C. Using (lC.2p and (1C.3P while keeping only leading terms, Eq. 
(IC.lj) becomes: 

^ = % M/* -*;)] + -^p (C.4) 

where t = t/L 2 and g has been replace by e _/3 / L ~ 1 — 4. The first term of RHS of Eq. 
(1C.4|) accounts for the drive which favors an ordered phase, whereas the second term 
represents the diffusion which favors a homogeneous phase. In the weak asymmetry 
limit, the two terms are comparable in magnitude and thus compete one another. 

These hydrodynamic equations are in fact coupled Burgers equations, whose 
stationary solution obeys 
dp A 



dx 
dp B 
dx 
dpc 



~P [pA (pB - Pc)} 

-f3 [ PB (p c - pa)} (C.5) 

-P [PC {pA - Pb)] , 



dx 

obtained by setting the LHS of Eq. (1C.4I) to zero and integrating over x. The absence 
of integration constant is due to the fact that there are no steady-state currents in the 
case of equal densities. 



Eqs. (1C.5P have been solved by Ayyer et al. 36j for the ABC model on an interval. 
For equal densities their derivation applies for periodic boundary condition as well. 
Multiplying the three equations by PbPc-, PaPc an d PaPb, respectively, and summing 
the resulting equations, yields -^(paPbPc) — 0, and consequently 

Pa{x)pb{x)pc{x) = K, (C.6) 



where K > is a constant. Equation ( 1C6|) in conjunction with Pa + Pb + Pc — 1 



decouples Eqs. (1C.5I) . yielding an equation for pa '■ 
dpA 



<)r ±PsJp A (l-p A ) 2 -±Kp A . (C.7) 

and similarly for pb and pc- In terms of the rescaled variables t = 2j3x and y(t) = pa{x) 
Eq. (1C.7I) may be written as 

2 

+ U K [y{t)\ = 0, (C.8) 



1 

2 

where 



dyjt) 
dt 



U K (y)= l -Ky- l -y 2 (l-yf. (C.9) 

This equation can be viewed as an equation of motion of a zero-energy particle with mass 
1 in a quartic potential. The four roots of Ux{y) = are denoted here as {0,a, b, c}. 
They are functions of K obeying 0<a<b<l<c. The physical trajectory, where 
< y{t) < 1 and Ux{y) < 0, is the one where the particle oscillates between a and b 
with a period of 

b 

T = 2 [ ; dy (CIO) 
J ^-2U K {y) 
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Based on the free energy of the ABC model, it can be shown that Eq. (10.81) has 
a unique steady-state solution given by T = 2/3, and many quasi-stationary solutions 
given by T = 2/3m where m > 1 is an integer. The trajectories are unique up to the 
choice of initial time for the motion of the particle which corresponds to the translation 
symmetry of the profile. In order to obtain an analytic expression for the integral in Eq. 
(IC.lOP it is convenient to rewrite it as an elliptic integral of the first kind of the form 

X 

r dz 

F(x,k)= , . (C.ll) 

V ; J v/(l - z 2 ) (1 - k 2 z 2 ) V ' 

This is done using a Mobius transformation that takes the roots of the potential, Ux{y), 
from {0, a, b, c} to {— 1, — 1/k, 1/k, 1}, which are the roots of the denominator of ( IC.lip . 
The transformation is given by 

» = /<V) = ^, (C12) 
a_y + 1 

where 

±ab + y/gb(c-b)(c- a) 
a± = aVc ' (CU3) 

and 

k = i±^£. (C.14) 
1 — a + a 

The parameters and k are functions of K through a, b and c. Let t(y) be the 

time it takes the particle to move from a to y. Using the transformation above it may 
be expressed as 

y 

t = 2> dy ' 



-2U K (y>) 

M 



"I w^m =>t[F(llk ' k)+F(mM ' (c ' 15) 



-1/k 



where 



K = —= :. (C.16) 

V (1 — a+a)(l — a+b)(l — a + c) 

The full period is given by setting f(y) = 1/k. The condition of T = 2(3 yields an 
equation that connects and K through k(K) and x(K): 

f3 = 2kF (l/k, k) . (C.17) 

With K(f3) known, we can proceed to express the profile by inverting Eq. (1C15[) . 
This is done using the Jacobi elliptic function, sn(x,k), defined by the equation 
F (sn (x, k) ,k) = x. The resulting profile is given up to translations of x as 

= l + m( 2g»/xt) 
a + — a„sn \2pxj x, k) 
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The profiles for the two other species are Pb{x) = p A {x — |) and pc{x) = Pa{x + |). 

The solution above of the standard ABC model with equal densities can easily be 
extended to the conserving model which includes vacancies. This is done by applying the 
mapping between the steady state of an L-size generalized ABC model to a 'condensed' 
iV-size system without vacancies but with an inverse temperature of (3r (see the first 
paragraph of Section H]). The mapping back to the L-size system requires the addition of 
vacancies into the lattice, done by multiplying the profile by r. Hence, the steady-state 
profile of the conserving model is given as 

1 + sn (20rx/ x, fc) 

Pa{x) = r — — - — — , (C.19) 

a + — a_sn (2prx/x, k) 

where x, k , a + and a_ are now functions of K((3r), and the two other profiles are again 
given as p B (x) = p A (x - §) and p c {x) = p A (x + \). 

We now consider the steady-state profile of the nonconserving model. The 
corresponding mean-field are obtained in a similar way by including the nonconserving 
process in Eq. (jSj), yielding 

^ = h+pL 2 I 2 , (C.20) 

where Ii = (3-^ [p A (pb — Pc)] + -q§t represents the drive and diffusion, and I 2 = 
Po ~ e ~ PaPbPc is corresponds to the evaporation and deposition processes. The 
equations for B and C are again given by cyclic permutations of this equation over 
A, B and C. Detailed balance with respect to the conserving (JTJ) and nonconserving 
dHJ) processes implies that in the steady state I\ = and I 2 = independently. Setting 
I\ = yields the same stationary equation and thus the same profile as in the conserving 
model. Since this profile obeys Pa(x)pb(x)pc{x) = Kr 3 and po(x) = 1 — r, 7 2 can indeed 
be set to zero for all x, yielding a relation between r and p : 



P = ~^n 



(C.21) 



This seemingly accidental coincidence is due to the choice of a nonconserving process 
that maintains detailed balance. 

Because both the conserving and nonconserving models have the same stationary 
profiles, we can use Eq. (1C.21|) to compute the chemical potential in the conserving 
model. For certain values of r we find that this definition yields a negative 
compressibility, ^ < 0. This profile is thus unstable in the nonconserving model, 
giving rise to inequivalence of ensembles. 
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